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We present a decomposition scheme based on Lie- Trotter-Suzuki product formulae to represent 
an ordered operator exponential as a product of ordinary operator exponentials. We provide a 
rigorous proof that does not use a time-displacement superoperator, and can be applied to non- 
analytic functions. Our proof provides explicit bounds on the error and includes cases where the 
functions are not infinitely differentiable. We show that Lie- Trotter-Suzuki product formulae can 
still be used for functions that are not infinitely differentiable, but that arbitrary order scaling may 
00 ■ not be achieved. 
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Q , I. INTRODUCTION 



Decompositions of operator exponentials are widely used to approximate operator exponentials that arise in both 
physics and applied mathematics. These approximations are used because it is often difficult to exponentiate an 
operator directly, even if the operator is a sum of a sequence of operators that can be easily exponentiated individually. 
The goal in a decomposition method is to approximate an exponential of a sum of operators as a product of operator 



exponentials. Particular examples of decompositions include the Trotter formula, and the related Baker-Campbell- 
Hausdorff formula as well as Suzuki decompositions [j], These approximations have found use in many fields 
including quantum Monte-Carlo calculations quantum computing Q and classical dynamics [ij to name a few. 
g The problem is to solve for the operator U given by 

d x U(X, n) = H(X)U(X, n), with U(jt, n) = 1, (1) 

■ where H is a linear operator and A and /i are real numbers. In general, computing U can be difficult regardless of 
whether H is dependent on A. The case of A-dependence makes the problem significantly more complicated and is 

. the focus of this work, whereas the case of A-independent H has been well studied. The Lie- Trotter formula gives 
a simple solution, whereas more efficient higher-order solutions are given by the Lie- Trotter-Suzuki (LTS) product 
O; formulae [li]. 

CNJ ■ A direct approach to solving U in the A-independent case is first to diagonalize H, then solve the differential 
equation (|TJ) by direct exponentiation. The scenario we consider is that this is not possible, and instead we are given 
a set of to operators {Hj : j = 1, . . . , to}, where it is possible to exponentiate of each of these operators, and 

H = J2 H r (2) 

X ■ 

The evolution operator can then be approximated by a product of exponentials of Hj for some sequence of {ji} and 
intervals {AAi}, 

N 

U(\,fi) « exp(F ijv AA JV )exp(ff ijv _ 1 AA^_i) • • • exp(H h AA X ) = ~[[e H » AX '. (3) 

t=i 

The goal is to make an intractable calculation of U tractable by approximating U as a finite-length product of efficiently 
calculated exponentials. The complexity of the calculation can then be quantified by the number of exponentials N, 
and the scaling of N in terms of the parameter difference AA = A — fi. 

The case of A-dependence makes the problem harder because, rather than the usual operator exponential of H, 
the solution is an ordered exponential. Diagonalization techniques are not directly applicable to solving ordered 
exponentials, and methods such as using the ordered product of exponentials are needed. That is, U is approximated 
by an expression of the form 



c/(A lA1 )«n eff(Ai)AAi - w 
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We consider the case where H may be A-dependent and where H is a sum as in ([2]). One approach would be to 
replace Eq. ([3]) with a product formula of ordered exponentials. There would still remain the problem of evaluating 
the ordered exponentials, which would require an approach such as (j4|). A simpler approach is to use j3|), but choose 
appropriate values of A at which to evaluate the Hj i . The approximation is then 

N 

t/(A, M )«n e ^ (A ° AAi - (5) 
i=i 

Suzuki provides an efficient method for approximating ordered exponentials in this way 0]. The method given by 
Suzuki is in terms of a time-displacement operator, but is equivalent. For many applications it is desirable to be able 
to place upper bounds on the error that can be obtained. Suzuki derives an order scaling, but not an upper bound 
on the error. Berry et al. find an upper bound on the error, but only in the case without A-dependence, and for H 
antihermitian (corresponding to Hamiltonian evolution) [j| . 

Here wc prove upper bounds on the error in the general case where H can depend on A, and is not restricted to be 
antihermitian. Whereas Suzuki uses a time-displacement operator, we provide a proof entirely without the use of this 
operator. The time-displacement operator is problematic, because it is unclear how it acts for A-dependence that is 
non-analytic. We find that, provided all derivatives of the Hj(X) exist, Suzuki's result holds. If there are derivatives 
that do not exist, then Suzuki's result does not necessarily hold; we demonstrate this via a counterexample. We solve 
the case where derivatives may not exist, and find that Suzuki's approach can still give better scaling of N with AA, 
although the scaling that can be obtained is limited by how many times the Hj(X) are differentiable. 

In Sec.|ll]we give the background for Trotter product formulae in detail. In Sec. lIIII we review Suzuki's decomposition 
methods, and provide our form of Suzuki's recursive method. In Sec. IIVI we introduce our terminology and present 
our main result. Then we rigorously prove the scaling of the error in Sec. [Vj and place an upper bound on the error 
in Sec. IVH We then use the error bounds in Sec. I VIII to find the appropriate order of the integrator to use. 



II. TROTTER FORMULAE 



Typically there are two different scenarios that may be considered. First, one may consider a short interval AA; 
the goal is then to obtain error that decreases rapidly as AA — * 0. Alternatively the interval AA may be long, and 
the goal is to obtain an approximation to within a certain error with as few exponentials as possible. For example, 
given A-independent operators A and B, it holds that 

e AX(A+B) = e AXA e AXB + (6) 

This gives an accurate approximation for small AA. For large AA, we may use Eq. © to derive the Trotter formula 

e AX(A+B) = (e AAA/„ e AAS/ n) „ + ( A A 2 /n). (7) 

This order of error is obtained because the error for interval AX/n is 0((AA/n) 2 ). Taking the power of n then gives n 
times this error if the norm of exp[(yl + B)AX] is at most one for any AA > 0, resulting in the error shown in Eq. 0. 
To obtain a given error e, the value of n must then scale as (9(AA 2 /e). The goal is to make the value of n needed to 
achieve a given accuracy as small as possible (n is proportional to the total number of exponentials). 

More generally, for a sum of an arbitrary number of operators Hj , similar formulae give the same scaling. To obtain 
better scaling, one can use a different product of exponentials. The Lie- Trotter-Suzuki product formulae [2[ replace 
the product for short AA with another that gives error scaling as 0(AA p+1 ). 

It can be seen that splitting large AA into n intervals as in Eq. yields an error scaling as 0(AX p+1 /n p ) if the 
norm of U is at most one for any A. It may at first appear that this gives worse results for large AA due to the higher 
power. In fact, there is an advantage due to the fact that a higher power of n is obtained. The value of n required to 
achieve a given error then scales as 0(AA 1+1 /p /e 1/p ). Therefore, for large AA, increasing p gives scaling of N that is 
close to linear in A A. 

Similar considerations hold for the case of ordered exponentials (i.e. with A-dependence). Huyghebaert and De 
Raedt showed how to generalize the Trotter formula to apply to ordered operator exponentials Q . Their formula has 
a decomposition error that is 0(AA 2 ), but requires that the integrals of A(u) and B(u) are known. Subsequently 
Suzuki developed a method to achieve error that scales as 0(AA p+1 ) for some ordered exponentials Q, and does not 
require the integrals of A and B to be known. We find that, in contrast to the A-independent case, it is not necessarily 
possible to obtain scaling as 0(AA p+1 ) for arbitrarily large p. It is possible if derivatives of all orders exist. If there 
are higher-order derivatives that do not exist, then it is still possible to use Suzuki's method to obtain error scaling 
as 0(AA p+1 ) for some values of p, but the maximum value of p for which this scaling can be proven depends on what 
orders of derivatives exist. 
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III. SUZUKI DECOMPOSITIONS 



In this section we explain Suzuki decompositions in more detail. In general, decompositions are of the form, as in 



JV 



Ufa + AA, ,x) = Y[ e ff * (A * )AA ' . (8) 
»=i 

Here we write the final parameter A as [i + A A, to emphasize the dependence on AA (= A — /x). There are many 
different types of decompositions, but the type that we focus on in this paper is symmetric decompositions because 
all Suzuki decompositions are symmetric. 

Definition 1. The operator U (jj, + AA, fx) is asymmetric decomposition of the operator U \fi + AA, /i) if U (/x + AA, /i) 
is a decomposition of Ufa + AA, xx) and Ufa + AA, /i) = [U fa, /x + AA)] -1 . 

An important method for generating symmetric decompositions is due to Suzuki 0, Q, which we call Suzuki's 
recursive method due to its similarity to the method presented by Suzuki in [2J . Furthermore we call any decomposition 
formula that is found using this method a Suzuki decomposition. 

Suzuki's recursive method takes a symmetric decomposition formula U v fa + AA,/i), that approximates an ordered 
operator exponential Ufa + AA,/x) with an approximation error that is at most proportional to AA 2p+1 as input, and 
outputs a symmetric approximation formula U p +\ fa + AA, /x) with an error that is often proportional to AA 2p+3 . The 
approximation U p +\fa + AA, /x) is found using the following recursion relations, 

Up+idi + AA, /x) ee U p (fi + AX,fi+ [1 - 8p]AX)U p (fj, + [1 - a p ]A\,n+[l - 2s p ]AX) 

x U p fa + [1 - 2s p ]AA, n + 2s p AX)U p fa + 2s p AX, /j + s p AA)C/ p (/i + s p AA, fjt), (9) 

with s p = (4-4 1 /( 2 p+ 1 ') _1 . 

Suzuki's recursive method does not actually approximate U(fi + AA, fi) but rather it builds a higher order approx- 
imation formula out of a lower order one. Therefore this method can only be used to approximate U(fi + AA, fi) if it 
is seeded with an appropriate initial approximation. A convenient approximation formula based on Suzuki's recursive 
method is the fc th order Lie- Trotter-Suzuki product formula which is defined as follows. 

Definition 2. The k th order Lie- Trotter- Suzuki product formula for the operator H(u) = Y)j—i Hj (u) and the interval 
\p, /j + AA] is defined to be Uk(fJ- + AA, fj,), which is found by using 



Ux{n + AA,ju) = ( Y[exp(H j (fj, + AA/2)AA/2) ] [ JJ exp(if,(/i + AA/2)AA/2) 

\J=1 / \3=m j 



(10) 



as an initial approximation and by applying Suzuki 's recursive method to it k — 1 times. 

Based on Suzuki's analysis 0, H[ Uk should have approximation error that is proportional to AA 2fc+1 . Hence if AA 
is sufficiently small, then the formula should be highly accurate. One might think that it would be advantageous 
to increase k without limit, in order to obtain increasingly accurate approximation formulae. This is not the case, 
because the number of terms in the formula increases exponentially with k. The best value of k to use can be expected 
to depend on the desired accuracy, as well as a range of other parameters [f| . 



IV. SUFFICIENCY CRITERION FOR DECOMPOSITION 



Suzuki's recursive method is a powerful technique for generating high-order decomposition formulae for ordered 
operator exponentials. The k th order Lie- Trotter-Suzuki product formula in particular seems to be well suited for 
approximating ordered operator exponentials that appear in quantum mechanics and in other fields; furthermore it 
appears that these formulae should be applicable to approximating the ordered exponentials of any finite dimen- 
sional operator H. However it turns out that Suzuki's recursive method does not always generate a higher order 
decomposition formula from a lower order one. 

We show this using the example of the operator H2{u) = u 3 sin(l/u)l. For this operator the second order Lie- 
Trotter-Suzuki product formula is not an approximation whose error as measured by the 2-norm is 0(AA 5 ). In Figure 
Q]we see that the error is proportional to A A 4 for the operator H a (u), rather than the A A 5 scaling that we expect 
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AA, A A, 

(a) (b) 

FIG. 1: This is a plot of C = \\U(AX,0) - U 2 (A\, 0)|| 2 /AA 5 for H a = w 3 sin(l/u)l in (a) and H b = cos(ii)H in (b). The error 
in (a) is proportional to AA 4 as opposed to the O(AA) 5 scaling predicted for that Suzuki decomposition. The error in (b) is 
proportional to AA 5 as expected for that Suzuki decomposition. 



and observe for the analytic operator Hb(u) = cos(u). This shows that the second order Lie- Trotter-Suzuki product 
formula is not as accurate as may be expected for some non-analytic operators. 

Our analysis will show that this discrepancy arises from the fact that H2(u) — u 3 sin(l/u) is not smooth enough 
for the second order Lie- Trotter-Suzuki formula to have an error which is (3(AA 5 ). In the subsequent discussion we 
will need to classify the smoothness of the operators that arise in decompositions. We use the smoothness criteria 
2/c-smooth and A-2/c-smooth, which we define below. 

Definition 3. The set of operators {Hj : j = 1, . . . , to} is P -smooth on the interval \p, A] if for each Hj the quantity 
\\d^ Hj(u)\\ is finite on the interval [pi, A]. 

Here, and throughout this paper, we define || • || to be the 2-norm. Also if {Hj} is P-smooth for every positive 
integer P, we call {Hj} oo-smooth. 

This condition is not precise enough for all of our purposes. For our error bounds we need to introduce the more 
precise condition of A-P-smoothness. This condition is useful because it guarantees that if the set {Hj} is A-P-smooth 
and p < P then < A p . This property allows us to write our error bounds in a form that does not contain 

any of the derivatives of H individually, but rather in terms of A which upper bounds the magnitude of any of these 
derivatives. We formally define this condition below. 

Definition 4. The set of operators {Hj : j = I, . . . , m} is A-P-smooth on the interval [/i, A] if {Hj} is P-smooth 
and A > sup p=0)lv .. iP ^sup„ e[ ^ A] (j2"U 11-^ Mil) ( J ■ 

For example if {H(u)} — {sin(2w)l}, where 1 is the identity operator, then using Definition g] {H(u)} is 2 2 / 3 - 
2-smooth on the interval [0, n] because the largest value \\H(u)^ \\ takes is 2 2 / 3 , for p — 0,1,2. It is also 

2-2-smooth because 2 2 / 3 < 2, furthermore since 

||#( u )(p)|ji/(p+i) < 2 for all positive inte gers p then {H(u)} is also 

2-oo-smooth. 

Using this measure of smoothness we can then state the following theorem, which is also the main theorem in this 
paper. 

Theorem 1. // the set {Hj} is A-2k-smooth on the interval [/i, pL + AA], and e < (9/10)(5/3) fc AAA, and e < 1 and 
max I>!/ ||(7(x,?/)|| < 1, then a decomposition U(p + AA, p) can be constructed such that \\U — U\\ < e and the number 
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of operator exponentials present in U , N, satisfies 

N < 2m5 k ~ 1 



l/2k 



(11) 



We prove Theorem [T] in several steps, the details of which are spread over Sees. IVland I VII In Sec. [V] we construct 
the Taylor series for an ordered operator exponential, and use this series to prove that the Lie- Trotter-Suzuki product 
formula can generate an approximation whose error is <3(AA 2fc+1 ), if {Hj} is 2fc-smooth on the interval [,u, n + AX]. 
In Sec. I VII we use the order estimates in Sec. |V]to obtain upper bounds on the error. The result of Theorem [T] then 
follows by counting the number of exponentials needed to make the error bound less than e. Finally in Sec. I VIII we 
show that if k is chosen appropriately, then N scales almost optimally with AA if there exists a value of A such that 
{Hj} is A-oo-smooth on [fi, /j, + AX] for every AA > 0. 



V. DECOMPOSING ORDERED EXPONENTIALS 



In this section we present a new derivation of Suzuki's recursive method. Our derivation has the advantage that it 
can be rigorously proven that if {Hj} is 2fc-smooth, where H (u) = 53j=i Hj(u), then the k th order Lie- Trotter-Suzuki 
product formula will have an error of 0(AA 2fc+1 ). We show this in three steps. We first give an expression for the 
Taylor series expansion of a ordered exponential U(/j, + AX,fi). Then using this expression for the Taylor series, we 
show in Theorem [2] that Suzuki's recursive method can be used to generate approximations to U((i + AA, /x) that 
invoke an error that is 0(AA 2fc+1 ) if {Hj} is 2fc-smooth. Finally we show in Corollary Q] that if {Hj} is 2fc-smooth 
then the k th order Lie- Trotter-Suzuki product formula has an error that is at most proportional to AA 2fe+1 . 

It is convenient to expand U in a Taylor series of the form 



U(p + AA, /i) = l+T 1 ( t i)AX 



T 2 (y)AX 2 
2! 



(12) 



If H is not analytic, then this Taylor series must be truncated, and the error can be bounded by the following lemma. 
Lemma 1. For all H : R — > C NxN , that are for P £ No, P times differ entiable on the interval [fi, fi + AX] C R, then 

p 



U(fi + AA, fj,) 



(AX) P T P 



p=0 



pi 



< 



{ ug[m,m+aa] \\Tp+i(u)U(u, n)\\AX p+1 



(13) 



where T p (u) is defined by the recursion relation, T p+ i(u) = T p (u)H(u) +dtT p (u), with Tq = 1 chosen to be the initial 
condition. 



Proof. 

We first show, for the positive integer £ < P + 1 that if s G [fi,fj,+ A A] then 



ds? 



U(8,n)=T t (s)U(8,ii). 



This equation can be validated by using induction on I. The base case follows by setting I = in |T 
demonstrate the induction step by noting that if (fT4)) is true for I < P then 



d 



+i 



d 



dse+T U(s,Li) = —T t {s)U{8,n) 
d 



ds 



Us) 



U(s, f i)+T e (s) 



ds 



U(a,n) 



(14) 
We then 



(15) 



Since Tg(s) contains derivatives of H(s) up to order I — 1, it follows that d s Tp(s) contains derivatives up to order 
£. Then since H(s) is P times differentiable, d s Tg(s) exists if I < P, which implies that 9f +1 C/(s,/i) exists. We 
then use the differential equation in (JTJ) to evaluate the derivative of U(fi + AX, fi) in (fT5]) and use the fact that 
T l+ i = HT e + d s T £ to find that 



gi+l 

■q^U{s,^i) = T l+ i{s)U(s,fi). 



(16) 
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This demonstrates the induction step in our proof of (fT4| . Since we have already shown that ([14]) is valid for To, it 
is also true for all Ti if I < P + 1 by induction on I. 

We then use (fl4|) and Taylor's Theorem to conclude that 



p=0 



(AX) p T p (p) 



We rearrange this result and find that 



U(fx + AA, fi) 



E 

p=0 



(AX) P T P 



pi 



< 



aa (AA-s) p 
Tp+i(fi + s)U(n + s,/i) As. 



max„ e [ M , M+ AA] \\Tp+i(u)U(u, fi)\\AX p+1 



(P + l)! 



(17) 

(18) 
□ 



Lemma Q] provides a convenient expression for the terms in the Taylor series of U(fi+ AA, fi), and it also estimates the 
error invoked by truncating the series at order P for any P £ N. The following theorem uses this Lemma to show that 
Suzuki's recursive method will produce a higher order approximation from a lower order symmetric approximation, 
if {Hj} is sufficiently smooth on [fi, fi + AA] and if H is the sum of all of the elements in the set {Hj}. 

Theorem 2. If H = X^Jli Hj where the set {Hj} is for a fixed p, 2(p + l)-smooth on the interval [fi, fi + AA] and 
U p (/i + AA, fi) is a symmetric approximation formula such that \\U p (fi + AA, fi) — U(fi + AA, /i)|| £ 0(AA 2p+1 ), and 
U p +\{fi + AA, fi) is found by applying Suzuki's recursive method on U p (pL + AA, /i) 7 then 



\\U(fi + A\,ti)-U p+1 (n + A\,(j,)\\ e 0(AX 2 p+ 3 ). 



(19) 



Proof. 

In this proof we compare the Taylor series of U to that of U p and show that by choosing s p appropriately will cause 
both the terms proportional to AA 2p+2 and AA 2p+3 to vanish. 

By expanding the recursive formula in Lemma [T] we see that a Taylor polynomial can be constructed for U whose 
difference from U is 0(AA 2p+3 ) because {Hj} is 2(p + l)-smooth on + AA]. A similar polynomial can be 
constructed for U p by Taylor expanding each Hj that appears in the exponentials in U p , and then expanding each 
of these exponentials. Then because {Hj} is 2(p + l)-smooth, Taylor's Theorem implies that this polynomial can be 
constructed such that the difference between it and U p is 0(AA 2p+3 ). Therefore since \\U — U p \\ S (9(AA 2p+1 ) there 
exist operators C and E that are independent of AA, such that 



U p (pl + AA, n) - U(n + AA, ^) = CV)AA 2p+1 + £(^)AA 2p+2 + 0(AA 2p 



+3\ 



(20) 



We then use the above equation to write U p +\ as 



U{ f i + AA, fi + [1 - s p ]AA) + C(/i + [1 - s p ] AA)(s p AA) 2p+1 + • ■ • ) x 



x[U(p + s p AX, n) + C(p)(s p AX) 



2p+l 



(21) 



Since {Hj} is 2{p + l)-smooth, and since H — Y^JLi Hj' ^ follows that U p +i is differentiable 2(p + 1) times. Then 
since U is differentiable 2(p + 1) times it follows from Taylor's theorem that C is differentiable, and hence we can 
Taylor expand each C in this formula in powers of AA to lowest order. By doing so and by defining E(fi) to be the 
sum of all the terms that are proportional to AA 2p+2 in this expansion we find that 



U p+1 (fi + AX,/j) 



U(fi + AA, pi) + [4s 2p+1 + [1 - 4s p ] 2p+1 ]C(/i)AA 2p+1 + E{p)AX 2p+2 + 0{AX 2p+3 ). 



(22) 



Then we see that if s p — (4 — 4 1 /( 2 f > + 1 )) 1 j then the terms of order 2p + 1 in the above equation vanish. Hence the 
error invoked using U p +i instead of U is 0(AA 2p+2 ) with this choice of s p . 

Next we show that E(/j,) = using reasoning that is similar to that used by Suzuki in his proof of his recursive 
method for the case where H is a constant operator [l|, Q. Because U p +\ is symmetric it follows from Definition [1] 
that 

1 = U p+ i(n, fi + AA)l/p + i(/i + AA, fi) 
= (u(ji, m + AA) + E(fi + AA)AA 2p + 2 ) (u(n + AA, /i) + J B( M )AA 2p+2 ) + 0(AA) 2p+3 

= 1+ \u(n,fi + AX)E(fi)+E(n + AX)U(fi + AX,fi)] AX 2p+2 + 0(AX) 2p+3 . (23) 
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This equation is only valid if 



U(fx, [i + AX)E(n) + E(ji + AX)U(fi + AX, n) e O(AA) 



We then show that E is zero by taking the limit of the above equation as AA approaches zero. But we need to 
ensure that E is continuous to evaluate this limit. The operator E consists of products of derivatives of elements from 
the set {Hj}, and these derivatives are of order at most 2p+ 1. Then since each Hj is diffcrcntiablc 2p + 2 times E 
is diffcrcntiablc, and hence it is continuous. Then using this fact it follows that 

A lim Q U(fx,fj, + AX)E(fj,)+E(n + AX)U(fj. + A\,n) =2E(n)=0. 

This implies that the norm of the difference between U and U p +\ is proportional to AA 2p+3 , which concludes our 
proof of Theorem [51 □ 

Now that we have proved that Suzuki's recursive method will generate a higher order decomposition formula from 
a lower order one if H is the sum of the elements from a sufficiently smooth set {Hj}, we now show that using the 
fc th order Lie- Trotter-Suzuki product formula invokes an error that is proportional to AA 2fc+1 if {Hj} is 2fc-smooth. 

Corollary 1. LetH{u) — Y^jLi Hj( u ) where the set {Hj} is Ik-smooth, on the interval [//, pt+ AA] and let U(fi+AX, /i) 
be the ordered operator exponential generated by H. Then if fi + AX) is the k th order Lie- Trotter- Suzuki formula, 
then 

||[/(/x + AA,/i)-[/ fe (M + AA,/i)|| e 0(AA 2fe+1 ). (24) 

Proof. 

Our proof of the corollary follows from an inductive argument on k. The validity of the base case can be verified by 
using Lemma[TJ More specifically, since {Hj} is 2fc-smooth on + AA] and since fc > 1, then H is at least three 
times differentiable on that interval. This means that we can use Lemma [T] to say that 

U(fi + AX,fi) = 1 + H(p)AX + [i? 2 (/i) + H'(n)]AX 2 /2 + 0(AA 3 ). (25) 

This expansion is also obtained by Taylor expanding exp(H(/i + AA/2)AA) to third order, so 

\\U(fi + AA, fi) — exp(H(fx + AA/2)AA)| e O(AA) 3 . (26) 

Since Ui(n + AX, /j,) is the Lie- Trotter formula for a constant H equal to H(p, + AA/2) it follows that, 

\\U 1 (fi,fi + AX) -exp(H(n + AA/2)AA)|| e O(AA) 3 . (27) 

It follows from the above equations and from the triangle inequality that the norm of the difference between U\ and 
U is at most proportional to AA 3 . 

Since we have shown that U±(fi + AA, /x) is a symmetric approximation formula whose error is 0(AA 3 ), it then 
follows from Theorem [2] and induction, that if {Hj} is 2fc-smooth then a symmetric approximation formula whose 
error is 0(AA 2fc+1 ) can be constructed from t/i(/x + AA, p) by applying Suzuki's recursive method to it k — 1 times. □ 

We have shown in this section that if {Hj} is 2&-smooth on the interval [fi, fi + A A] and if p < k then Suzuki's 
recursive method can be used to create a symmetric decomposition whose error is 0(AA) 2p+1 out of a symmetric 
decomposition whose error is 0(AA) 2p_1 . Then we have used this fact to show that the norm of the difference between 
U([j, + AX, fi) and the fc th order Lie- Trotter-Suzuki formula is (9(AA 2fc+1 ). In the following section we strengthen this 
result by providing an upper bound on the error invoked by using the fc th order Lie- Trotter-Suzuki product formula. 



VI. ERROR BOUNDS AND CONVERGENCE FOR DECOMPOSITION 

We showed in Sec. [V] that if the k th order Lie- Trotter-Suzuki product formula is used in the place of the ordered 
operator exponential of H, then an error is incurred that is at most proportional to AA 2fe+1 if H = Hj and 

the set of operators {Hj} is sufficiently smooth. We also showed that a sufficient condition for smoothness of the 
set {Hj} is a condition that we called 2fc-smooth, where this condition is defined is Definition [3] In this section we 
extend that result by finding upper bounds on the error invoked in using the Lie- Trotter-Suzuki product formula 
to approximate ordered operator exponentials if {Hj} is A-2fc-smooth. Unlike the previous section, here we assume 
that ma,x x>y \\U(x,y)\\ is at most one. This assumption is important because it ensures that our error bounds are 
not exponentially large. Our work can be made applicable to the case where this norm is greater than one by 
re-normalizing U. We discuss the implications of this in Appendix B. 
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We first provide in this section an upper bound on the error invoked in using the k th order Lie- Trotter-Suzuki 
product formula to approximate the ordered operator exponential Ufa + AA, /x) if AA is sufficiently short. We then 
use this result to upper bound the error if A A is not short. More specifically, we show that for every 1 < e > and 
AA > there exists an integer r such that 



Ufa + AA, p) - [J U k fa + qAX/r, p+(q- l)AA/r) 



9=1 



(28) 



if {Hj(u)} is 2fc-smooth on the interval + AA]. Finally by multiplying the number of exponentials in each k 



th 



order Lie- Trotter-Suzuki product formula by r, we find the number of exponentials used in the product in (|28|). We 
then use this result to prove Theorem Q] 

Our upper bound on the error invoked by using a single Uk to approximate the ordered operator exponential 
Ufa + AA, p) is given in Theorem [3j Before stating Theorem [3] we first define the following terms. Since the 
k th order Lie- Trotter-Suzuki product formula is a product of 2m5 fc ~ 1 exponentials, we can express this product as 



n 



2mn" 



exp(Hj c fa c )AX c ). We then use this expansion to define the following two useful quantities. 



Definition 5. We define q c .2k = %x and also define Qk = max c |<7 c ,2fc|- 

It can be shown that Qi = 1/2 and that if p > 1 then Q p = |1 — 4si| ■ • • |1 — 4s p _i|. We show in Appendix A 
that for any integer p that Q p < 2p/3 p , implying that Q p decreases exponentially with p. We use this definition of 
Qk in the following Theorem, that gives an upper bound on the difference between the ordered operator exponential 
Ufa + AA, p), and the fc th order Lie- Trotter-Suzuki product formula Ukfa + AA, p). 

Theorem 3. Let H(u) = ^ {Hj( u )} be A-2k- Suzuki- smooth on the interval \p,(J> + AA], and let 

max 1>9 ||l/(x,2/)|| < 1. Then if 2V2(5) fc ~ 1 Q fc AAA < 1/2, it follows that 

\\U(n + AA, n) - C/ fc (M + AA, < 2 [3(5) fe - 1 Q fe AAA] 2 " +1 , 
where Uk is given in Definition^ 

The proof of Theorem [3] requires us to first prove two Lemmas before we can conclude that the theorem is valid. 
We now introduce some notation to state these lemmas concisely. Since we have assumed that {Hj} is 2fc-smooth, 
Theorem [2] implies that the difference between U(fi + AA, p) and [/&(/! + AA, p) is 0(AA) 2fc+1 . Then using this fact, 
we know that we only need to compare the terms of 0(AA 2fc+1 ) to bound the difference between U and Uk- We 
introduce the following notation to denote only those terms that do not necessarily cancel. 

Definition 6. If the operator A{ AA) can be written as A(AX) = Y^,p=o A,AA P + i?(AA) where the norm of i?(AA) is 
0(AA) 2fc+1 , then we define R 2fc [A(AA)] to be the norm ofR(AX). 

This definition simply means that R.2k is the error term for a Taylor expansion to order 2k. Then using this 
definition, it follows from the triangle inequality that the norm of the difference between U and Uk is at most 

R 2fc [Ufa + AA, /i)] + R 2fc [Ukfa + AA, p)} . (29) 

Our proof of Theorem [3] then follows from (f2T)|) and upper bounds that we place on H 2 k[U fa + AX, p)} and R 2 fc [£4 (/i + 
AA, p)]. Our bound on R2k[Ufa + AA, p)] follows directly from LemmaQ] but the bound on ~K 2 k\Ukfa + AA, p)] does 
not. We will provide the latter upper bound in Lemma [3j but first we provide Definition [7] and Lemma [2] 

Definition 7. Let k be an integer and let H = Y^JLi then Ukfa + AA, AA) can be written as a product of the form 



n 



exp(Hj c (p c )AX c ). We then define X p for p < 2k to be 



and for p — 2k we define X 2 k to be 



X,, 



2m5 

E 

c=l 



fac - M) P 

AXp 



\Qe,2k\ 



(30) 



X 



2k 



2m5 K ~ 

E 

c=l 



max 
Te[(i,/i+AA] 



Hf k \r) 



fac - M) 2 fc 

AA 2fc 



\Oc,2k\- 



(31) 



Here the quantity q c ^2k is given in Definition^ 
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Then using this definition our lemma can be expressed as follows. 

Lemma 2. Let H(u) = Y) ■_ 1 Hj{u) and let the set {Hj} be Ik-smooth on the interval [fi,fi-\- AX], then the norm of 
the difference between Uj~ (ft + AA, ft) and its Taylor series in powers of AX truncated at order 2k, is bounded above by 



R 



2/.- 



2 k 



-p ( E §^ p+1 



\P=0 



pi 



(32) 



Proof. We begin our proof of Lemma [2] by writing Uk as a product of 2m5 k 1 exponentials and use Taylor's theorem 
to write Uk as 



2m5 fc ~ 1 

n exp 



c=l 



'2fc-l tt(p) 



, p=0 



(2k), s(Mc-s) 



2fc-l 



(2fc- 1)! 



-ds A A c 



We introduce the terms w c = (/i c — fi)/ AX and = an d use them to write Uk(p- + AX, fi) as 



(33) 



2m', 



n exp 



v p=0 



p! 



gfV + xAX) {Vc {2k ^ (AA) 2fc ^ j |^|AA 



(34) 



We now prove the lemma by placing an upper bound on H 2 k[Uk(p + AA, fi)] by expanding this equation in powers 
of AA, while retaining only those terms of order 2k + 1 and higher. As mentioned previously, the lower order terms 
are irrelevant since Theorem [2] guarantees that they cancel. 

By expanding the exponentials in (|34|) . taking the norm, using the triangle inequality, upper bounding each of the 
norms present in the expansion, and collecting terms again, we find that an upper bound on R,2fe [£/&(// + AX, fi)] is 



R 2 k exp 



2»i5 



, »-i ||i/, M («) 

p=0 ^ 



u?AA p 



max 

TG[At,At+AA] 



4 2fc) (-) 



(^ c AA) 2 fc 
(2fc)! 



|9cfc|AA 



(35) 



This equation can be simplified by substituting the constants X p into it. These constants are introduced in Definition 
[7J After this substitution our upper bound becomes 



R 



2/,' 



2k 



-p(Et AAP+1 



(36) 



□ 



We use Lemma [2] to provide an upper bound on the sum of the norm of all terms in the Taylor expansion of 
C/fc (/^ + A A, fi) which are of order 2k + 1 or higher. This bound is given in the following lemma. 



Lemma 3. Let H(u) = J2T=i Hj{u) where the set {Hj} is A-2k- 



smooth on 



the interval [fi, [i + AA] and let 



2v / 2(5)' c 1 QfcAAA < \. Then the norm of the difference between + AA, \i) and its Taylor series in AX truncated 
at order 2k, is upper bounded by 



H 2 k[U k (p + AA, ft)] < 2 f2v / 2(5) fc - 1 Q fc AAA) 



2fe+l 



(37) 



Proof. 

To simplify the following discussion we introduce T 2 k, defined by 



^2k 



ji/(p+i]. 



p=0,...,2fc 

We first find an upper bound on T 2 k ■ Now, by Definition [Jj 

2m5 fc ~ 1 



x n < 



E 1^1 



tic- H 

AX 



max 
t£[ai,m+aa] 



(38) 



(39) 
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In Appendix A, we show that q c . k < Qk < ft for all c, where Q k is an upper bound on the q Ctk given in Definition [5] 
For a 2fc th order Lie- Trotter-Suzuki product formula, ^"a^ — 1 ■ P m S'g m § these two bounds into the above inequality 
and using the fact that each element of {Hj} occurs 2(5 fe ~ 1 ) times in U k yields, 



X p <2(5) fe - 1 Q fc V max 



t— f tS[m,M+AA] 



tf, (p) (-) 



(40) 



Since we assume that {if,} is A-2fc-smooth, then \\H { f ] < A?p +1 , so Xp /(p+1) < (2{5) k - 1 Q k ) 1/(p+1) A. We 

also show in Appendix A that Q k > ^^r, implying that 2(5) k ~ 1 Qk > 1, and thus that Xp /(p+1) < 2(5) fc - 1 Q fc A. 
Since this upper bound holds for all < p < 2k, we conclude that 



r 2fe < 2(5)*- x Q fc A 



We begin the main inequality in Lemma [3j by expanding R 2 fc exp ^X)p=o AA P+1 



(41) 

in powers of AA and using 



pS is ^p+i — ana - -Xp < 1 . Writing the resulting expansion as an exponential we find that 



R 



2/? 



2 k 



A. 



Vp=0 



< R 



2 A- 



(V2F 2fe AA)P+ 1> 

cx p 2^ ■ 

\p=0 



p + 1 



It then follows that the right hand side of the above expression is upper bounded by 

(V2F 2fe AA)P+ lS 



R 



2 k 



exp 



p + 1 



Using the Taylor expansion of ln(l — x), we rewrite this as 



(42) 



(43) 



R 



2/v 



cxp (- ln(l - v^FafcAA) 



R 



2/,- 



1 - \/2F 2fe AA 



53 (72r 2fc AAf. 

p=2k+l 



(44) 



Provided \/2F 2fe AA < ~, this is upper bounded by 2(V2F 2fe AA) 2 ' £+1 . Plugging inequality JHJ) into Eq. £44} then 
gives that 



R 



2/. 



2 A' 



\p=0 



< 2 (2V2(5) fe - 1 g fe AAA 



2fe+l 



(45) 



The lemma follows by applying Lemma O 



□ 



Now that we have proven Lemma[3]we have an upper bound on R 2 fc [[/&(/! + AA,/i)]. We now use this upper bound 
to prove Theorem [3l 



Proof of Theorem [31 

Our proof of Theorem [3] begins by recalling the fact that 

R 2fe [U(ji + AA, /*) - U k (/x + AA, /*)] < R 2fc [1% + AA, //)] + R 2fc [U k {p + AA, //)] 



(46) 



We then place an upper bound on R 2 /j [t/(/i + AA, (i)] using Lemma [T] Using the notation of Lemma [T] we write the 
Taylor series of C/(/i + AA, /z) as J2 P T p AX p /p\. We then use the assumption that \\U(fi + AA, /i) | is less than one, to 
show from Lemma [1] that R 2 ^ [U(fi + AA, /i)] is at most 



max Me[ ^ +AA] ||T 2fc+ i(u)||AA 2fc+1 
(2fc + l)! 



(47) 



Using the recursive relations in Lemma [T] it follows that T 2 fe+i can be written as a sum of (2k + 1)1 terms that are 
each products of H and its derivatives. Then since {Hj : j = 1, . . . , m} is A-2fc-smooth and H — Ylj=i Hj> ^ follows 
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from Definition H that \\H^(u)\\ < A? for all u in the interval [//, /i + AA]. It can then be verified that each term in 
T 2 k+i must have a norm that is less than A 2fc+1 . Therefore it follows that || T^fc+i || < (2fc + l)!A 2fc+1 and hence 



R 2fc [U(fi + AA,/x)] < (AAA) 



2fc+l 



(48) 



Using Eq. ([46]) and Lemma [3] we see that if 2V2(5) fc - 1 Q fe AAA < 1/2 then an upper bound on the sum of 
R 2 fe [U(jx + AA, n)} and R 2fc [J7 fe (/x + AA, /x)] is 



(AAA) 2fc+1 + 2[2V2(5) fc ~ 1 Q /£ AAA] 2 ' £+1 . 
We then replace this upper bound with the following simpler upper bound 

2[3(5) fc - 1 Q fe AAA] 2fc+1 . 
This is the claim in Theorem [3j and hence we have proven the theorem. 



(49) 

(50) 
□ 



The error bound in Theorem [3] is vital to our remaining work, because it provides us with an upper bound on 
the error invoked by approximating an ordered operator exponential by Uk(fJ- + AA, fi) if AA is short. We will now 
show a method to devise accurate approximations to the ordered operator exponential U(fj, + AA, fi) even if AA is not 
short. Our approach is similar to that used by Berry et al. in [f| and that used by Suzuki in we split the ordered 
exponential into a product of ordinary exponentials, each of which has a short duration. However to do so we need to 
present a method to relate the error invoked by using one Uk to the error invoked by using a product of them. This 
result is provided in the following lemma. 

Lemma 4. If \\A p — B p \\ < 6/P where 6 is a positive number less than 1/2 and \\A p \\ < 1 for every p S {1, 2, • ■ • , P} 

then the product \\ Y[p=i A p — Yip=i Bp\\ — 2<$- 

Proof. 

Our proof begins by assuming that there exists some integer q such that 



IK IK 

p—1 p—1 



< 



p 



(51) 



We then prove Lemma 0] by using induction on q. The proof of the base case follows from \\A p — B p \\ < S/P. We 
then begin to prove the induction step by noting that from \\A p — B p \\ < 5/ P there exists an operator C with norm 
at most one, such that B q+ i = A q+ \ + (5/ P)C. Then by making this substitution and using the triangle inequality 
it follows that 



9+1 9+1 

p—i p—i 



< 



A 



9+1 



n a p - n r p 

\p=i p=i 



(S/P) 



P =i 



(52) 



Then because \\A p \\ < 1 and ||£? p || < 1 + S/P it can be verified using our induction hypothesis that the left hand side 
of Equation (|52|) is bounded above by 



( g + 1)^(1 + -^ 
P 



(53) 



This proves our induction step, and so it follows that || Y[ P =i Ap — 11^=1 ^p\\ — ^(1 + &/P) P 1 by using induction on 
q until q = P. The proof of the Lemma then follows from the fact that if 5 < 1/2 then (1 + (5/P) p_1 < 2. □ 



Using Lemma 0] we can now place an upper bound on the error for decompositions with longer AA. 

Lemma 5. // H (u) = X)J=i Hj(u) is 2k- Suzuki- smooth on the interval [/i, /i + AA], and max I>!/ ||C^(a^, 
e < 3Q fc (5) fc - 1 AAA, where Q k is given in Definition and the positive integer r is greater than 

2(3Q fc (5) fc - 1 AAA) 1+1 / 2fc 

e l/2k 

then we obtain that 

r 

U(n + AA, /i) - IJ E4(/x + qAA/r, fi + (q - 1) AA/r) < e, 

9=1 

where Uk is the k th order Lie- Trotter- Suzuki product formula, which we introduced in Definition^ 



< 1 and 



(54) 



(55) 



12 



Proof. 

Using the bound e < 3Q fc (5)* _1 AAA, we find using Eq. ([SI]) that 3Q fc (5) fe " 1 AAA/r < 1/2. Hence we can use 
Theorem [3] to obtain, for each q = 1, . . . , r, 

\\U(ji + qAX/r, n + (q - l)AA/r) - U k (p + qAX/r, n+(q- l)AA/r)|| < 2(3Q fc (5)^ 1 AAA/r) 2fe+1 . (56) 

We then re-write this bound as, 

2(3Q fc (5) fc - 1 AAA/r) 2fe + 1 = g ^(S^AAA) 2 ** 1 _ (5?) 
Then from dHJ) we can see that, because r > 2(3Q fe (5) fe " 1 AAA) 1+1 / 2 ' £ /e 1 / 2 ' £ it follows that 

2(3Q,(5) fc - 1 AAA/r) 2fe+1 < ~ (^) . (58) 

Then since k > 1 it follows that 

||Z7(/* + gAA/r, M + (g - l)AA/r) -U k ( f M + qAX/r, f x + (g-l)AA/r)|| < ^. (59) 

2r 

Then since both e and max x>y ||{7(x,j/)|| are less than one, the result of this lemma follows from Lemma |4] □ 

Lemma [5] shows that if the maximum value of the norm of U is one, then a product of k th order Lie- Trotter-Suzuki 
formulae converges to U as r increases. Furthermore we can also use this result to prove Theorem[T]by using the value 
of r from this Lemma and multiplying it by the number of exponentials in each U k to find a bound on the number of 
exponentials that are needed to approximate U(n + AA,/i). This proof is presented below. 

Proof of Theorem^ It can be verified from the definition of f/fe(/x + AX, fi) in Theorem [5] that there are at most 
2m5 fc_1 exponentials in each U k and since at most r different Uk in are needed to approximate U within an error of 
e then if max x>y \\U(x,y)\\ and e are at most one, it follows from Lemma[5]that the number of exponentials used to 
decompose U(fx + AA, fi) is at most 



N < 2m5 fc ~V < 2m5 fe_1 



2(3(5) fc - 1 AQ fc AA) 1+1 /2fc 



,l/2fe 



(60) 



if € < 3Qfc5 fc - 1 AAA. We then use the upper bound from Appendix A, Q k < 2fc/3 fc and the fact that (2fc) 1 / 2fe < 1.5 
to show that, 



N < 2m5 



fc-i 



k / A A \ \ l/2fe' 



., t v , , 5V / AAX\ 
5fcAAA - 



(61) 



Equation l|6ip is only valid if e < 3Qfc5 fc_1 AAA, this bound can be simplified by using the lower bound on Qk in (|A2|) . 
After substituting this lower bound we then find that e < (9/10)(5/3)' c AAA also is sufficient to guarantee that (|6ip is 
valid, which proves our theorem. □ 

Theorem Q] provides an upper bound on the number of exponentials that are needed to decompose an ordered 
operator exponential using a product of fc th order Lie- Trotter-Suzuki product formula, while guaranteeing that the 
approximation error is at most e. In the following section we present a formula that provides a reasonable value of k, 
for a particular set of values for e, A and AA. Furthermore we show that if {Hj} is A-oo-smooth and if that formula 
for k is used, then the number of exponentials used scales near optimally with AA. 



VII. ALMOST LINEAR SCALING 



Reference [5[ shows that there exist operator exponentials that, when decomposed into a sequence of N exponentials, 
require that N scale at least linearly with A A for large AA. This implies that any decomposition method that does 
not use any special properties of the operator being exponentiated, will also require that scale at least linearly 
with AA. We now show that if there exists a A such that the set of operators {Hj} is A-oo-smooth on [fx, \i + A A] for 
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every AA > 0, then we can choose k such that N scales almost linearly in AA. Specifically, we show that N/AX is 
sub-polynomial in AA, i.e., that liniAA^oo A ^T +d = for all constants d > 0, provided that max I>9 y)|| < 1. 



It follows from Theorem Q] and property [a] < a + 1 for any positive a, that if AAA > 1 and e < 1 then 



< 3mfcAAA 



25 



/AAA\ 1/2fc 



We set 



so that (f ) 



fc = 



/l ( AAA A 

2 l0g25 / 3 l^J 



> 



' AAA 



l/2fco 



Then 

N 
AA 



< 3mfc n A 



25 x 2fe 



3™fcoAexp fco21n 



(62) 



(63) 



(64) 



which is sub-polynomial in AA, though not poly-logarithmic in AA. In conclusion, if H is A-oo-smooth and if we 
choose k = fc , then the number of exponentials needed to decompose a ordered operator exponential of H using the 
fc th order Lie- Trotter-Suzuki formula scales almost linearly in AA. 

This choice of fco will cause A~ to scale nearly linearly with AA if H is oo-smooth; however if {Hj} is only 2P-smooth 
for some positive integer P, then we do not expect this because Theorem [T] cannot be used if fco > P. Hence a 
reasonable choice of fco is 



fco = min < P, 



'1 /AAA 



(65) 



The choice of fco in (|65|) does not allow for near linear scaling of A?" with AA, but it does cause N to be proportional 
to AA 1+1 /( 2P ) in the limit of large AA, and causes N to have the same scaling with AA that a A-oo-smooth {Hj} 
would have if A A is sufficiently short. 

In this section we require that max I>!/ ||t/(x, y)|| < 1 for this near linear scaling result to hold, but if this inequality 
does not hold then U can be normalized to ensure that it does, so it may seem that this result is more general than 
we claim. However we note in Appendix B that the un-normalized error in the decomposition of the un-normalized 
U can vary exponentially with AA. As a result we can only guarantee that the value of N needed to ensure that 
\\U — U\\ < e can be chosen to scale near linearly with AA if max x>y \\U(x,y)\\ < 1. 



VIII. CONCLUSIONS 



We have presented in this paper a rigorous derivation of Suzuki's recursive method for generating higher order 
approximations to ordered operator exponentials from a lower order formula. We have also shown that if H(u) = 
y^,j—i Hj(u) and {Hj} is 2fc-smooth, a condition which we define in Defmition[3l then Suzuki's recursive method can 
be used to build approximation formulae for the ordered exponential of H while invoking an error that is at most 
proportional to A\ 2k+1 . Furthermore we have shown that the fc th order Lie- Trotter-Suzuki product formula has an 
error at most proportional to AA 2fc+1 if H is 2fc-smooth. 

We have also shown that if H is A-2fc-smooth, which is a condition that we give in Definition [U and if 
max. x>y \\U(x, y)\\ < 1 then the number of exponentials needed to approximate an ordered operator exponential 
of H using a product of fc th order Lie- Trotter-Suzuki formulae while invoking a total error of at most e is bounded 
above by 

Finally we have shown that if {Hj} is A-oo-smooth then for long simulations, the number of exponentials needed to 
approximate ordered operator exponentials of H with an error of at most e scales close to linearly with AA. Linear 
scaling has been shown to be optimal [5;], so in this sense our scheme is nearly optimal. 

An important extension of this work will be to provide upper bounds on the number of exponentials used to 
decompose an ordered operator exponential when the user chooses step size adaptively, rather than the constant sized 
steps that we use in our derivation. Choosing these steps adaptively could lead to substantial improvements in the 
performance of our decomposition for certain ordered exponentials. 



N < 2m5 k - 1 
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APPENDIX A: DERIVATION OF BOUNDS ON Q k 

When proving bounds on the error introduced by our decomposition of an ordered exponential in Section IVI1 we 
use the quantities Qk defined by 

Qk = ^max{s 1 ,|l-4s 1 |}max{s 2 ,|l-4s 2 |}---max{s fe _ 1 ,|l-4s fc _ 1 |} (Al) 

for k > 1 and Qi = -k, where Sk — zzzt7T2Ffti f° r all fe > 1. We now show that Qk decreases exponentially in fc, 



3 1 n 2k 

-—r < Qk< - 

2 3 fc ~ ^ 3 



k -<Q k < ^- (A2) 



The lower bound follows directly by noting that Sk > § for all k > 1. 

Set a — 21n(4), which is approximately 2.7726. Using that ~x > ln(l — x) for < x < 1, we then have that for 
k > 1, 

' M4) = TT7^^>ln(l-^ V (A3) 



2k + 1 w 2(2fc+ 1) - V 2(2fc + 1 
Taking exponentials on either side yields, 



4 -i/(2*+i) > l _ a (A4) 
2(2fc + l) V ; 



Multiplying by 4 and subtracting 1 on either side gives, 



and taking reciprocals then yields, 



1 1 k+h Ik 



l 



4Sfc 1_ ^k/(2k+i)_ 1 ^ 3fe + i-f " 3fc + i-l' (A6) 



Noting that Ask — 1 > s/c since Sk > i, and using that si < |, we conclude that 



for k > 2, and by inspection that the inequality Qk < also holds for k = 1 



12 1 k- i 2/c , , N 

Qk ~ 2 3 3^^ " F (A7) 

2fr 



APPENDIX B: NORMS LARGER THAN 1 



In this work we have restricted the norm of U(X,ii) to not exceed 1. This means that the eigenvalues of H(\) can 
have no positive real part. In the case where they do, then the analysis can be performed in the following way. Simply 
define the new operators 



H'(X) = H(X) - k(A)1, 
H' j (X) = H j (\)-( K (X)/m)t, 



(Bl) 
(B2) 
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for some k(X) such that the eigenvalues of H'(X) have no positive real part. Then the result we have given in Theorem 
Q]will hold for H' and {H'^} (provided we also define A in terms of these operators). The difference between H and 
H' simply corresponds to a normalization factor; i.e. 

t/(A,/x) = U'{\,p)eS> {x)dx . (B3) 
To approximate U(X, /i) by a series of exponentials, we can simply use the series to approximate U'(X, /x), which gives 

N 



i=l 

N 

= irJJexp( J ff ij (A i )AA i ), 



(B4) 



where K is a normalization correction 



K 



(B5) 



Thus the same series of exponentials can be used, except for a normalization factor. There is a difference in the 
final error that can be obtained, because 



^'(A, Ai )-[]exp(^(A 4 )AA l ) 



i=l 



< e 



implies that 



N 



f k(x)cLx 



(B6) 



(B7) 



It might be imagined that the relative error can be kept below e with similar scaling of N. That is, that e^ K ( x ) dx can 
replaced with ||C/(A, n)\\ in (|B7|) . Unfortunately, that is not the case. The reason is that, due to the submultiplicativity 

of the operator norm, ||J7(A, /x)|| can be much smaller than eU h[x> x . 

For example, consider the case where H is initially a z (the Pauli operator) over an interval AX/2, then is — a z 



over another interval AA/2. Then C/(A,/x) = 1, and has norm 1, but 



„/„ K(x)dx _ aa 



A small error in between the 



two intervals of length AA/2 can then yield a large relative error in the final result. For example, consider the error 
E = e lSa y . That yields a final result 



[/(A, it + AX/2)EU(jjl + AA/2, fx) 



cos S e AA sin S 
s AA sin S cos S 



(B8) 



The error in this result scales as e AA , despite the final norm being small for U(X, fj.). 

With the possibility that the norm of U(X, /i) exceeds 1, our approach need not give scaling for N that is close to 
linear in AA. In the lower bound on N in Theorem[T] the {l/e) 1 / 2k will be replaced with 



(B9) 



To prevent this term scaling exponentially in AA, one would need to take k proportional to AA. However, this would 
result in (25/3)* scaling exponentially in AA. As a result, it does not appear to be possible to obtain subexponential 
scaling if there is no bound on the norm of U(X, /z). 
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